Deterministic Generation of Entangled Photons in Superconducting Resonator Arrays 
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We present a scheme for the deterministic generation of entangled photon pairs in a supercon- 
ducting resonator array. The resonators form a Jaynes-Cummings Lattice via the coupling to super- 
conducting qubits and the Kerr-like nonlinearity arises due to the coupling. We show that entangled 
photons can be generated on-demand by applying spectroscopic techniques and exploiting the non- 
linearity and symmetry in the resonators. The scheme is robust against small parameter spreads due 
to fabrication errors. Our findings can be used as a key element for quantum information processing 
in superconducting quantum circuits. 
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Recent progresses in superconducting quantum circuits 
[H have enabled intensive study of the quantum be- 
havior of the microwave cavity modes in superconduct- 
ing resonators [2]. Strong coupling between the cavity 
modes and superconducting qubits has been experimen- 
tally demonstrated in a number of systems [3]. Quan- 
tum optical effects such as single photon manipulation 
and nonlinear spectrum have been observed [4]. Given 
their ultra-high quality factors, the superconducting res- 
onators can be a powerful platform for quantum informa- 
tion processing and quantum state engineering involving 
microwave photons. Recently, the Kerr nonlinearity was 
explored for the parametric amplification and squeezing 
of the resonator modes [5|,|6(- The generation of twin pho- 
tons at different frequencies was studied via the pumping 
of superconducting qubits [7]. Meanwhile, coupled cav- 
ity arrays (CCA) can also be realized in superconducting 
systems to study quantum many-body effects such as the 
Mott insulator-to-superfluid transition [8]. 

One essential element for quantum information process- 
ing using microwave photons is the entangled photon 
source that can generate entangled photon pairs in a 
deterministic way and to distribute the photons in the 
circuits [9(|. The entangled photons can be used to im- 
plement quantum teleportation and quantum cryptogra- 
phy and to test quantum mechanical principles such as 



the Bell inequalities [ICfl, UJJ. Various systems and ap- 



proaches were explored to generate entangled photons in 
the past, including spontaneous parametric down conver- 
sion (SPDC) in quantum optical systems |12j, pumping 
of bi-excitons and cavity-polaritons in semiconductors 
[13], and most recently, using the atom-photon interface 
in coupled cavities [14] . In those approaches, entangle- 
ment generation is either a stochastic process or a direct 
switching from the atomic to cavity states. 

In this work, we present a novel approach for the de- 
terministic generation of spatially-entangled or energy- 
entangled photons in superconducting resonators con- 
nected in a ring geometry. Microwave photons can tun- 
nel between adjacent resonators via capacitive coupling. 
Each resonator is coupled to a superconducting qubit to 
form the Jaynes-Cummings (JC) Lattice, which can give 



rise to the Kerr-like nonlinearity in the resonator modes 
(TBI [l6| . By combining the nonlinearity, circuit symme- 
try, and a spectroscopic technique, we show that entan- 
gled photons can be generated with high fidelity in the 
presence of resonator dissipation and qubit decoherence. 
The nonlinearity and the symmetry are exploited to pre- 
vent off-resonant transitions which would otherwise im- 
pair the entanglement generation. This idea can be fur- 
ther explored to generate novel quantum states involving 
multiple photons and can also be applied to semiconduc- 
tor micro-cavities. 

The proposed circuit includes four superconducting res- 
onators coupled in a ring geometry and each resonator 
couples to a superconducting qubit, as is shown in Fig.Q] 
The system forms a JC Lattice with the total Hamilto- 

nian H t = T, j ( H jc + H c ) ) which includes the JC model 
for each resonator-qubit system 



H 



^ = huoa^cij + (huo z /2)a jz + %(at<7j_ + Vj+CLj) (1) 
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and the tunneling between adjacent sites H t 
HJ(cLjaj+i + ctj^dj ). Here, Oj (aj) is the creation (an- 
nihilation) operator of the resonators, <jj z , <jj± are the 
Pauli matrices for the qubits, uo is the resonator fre- 
quency, hg is the coupling between the resonator and 
the qubit, and HJ is the tunneling matrix element be- 
tween the resonators. The eigenstates of the JC Model 
(i) 

Hjq are polariton states containing excitations in both 
the resonator and the qubit [2]. With a finite coupling 
and a small detuning A = uj z — c^o, the polariton spec- 
trum demonstrates Kerr-like nonlinearity similar to the 
onsite interaction in the Bose-Hubbard Model. In the 
thermodynamic limit, the JC Lattice exhibits the Mott 
insulator-to-superfluid phase transition d, [ItJ • 

The nonlinearity and the symmetry in the JC Lattice 
play an important role in the eigenenergy spectrum of the 
resonator array [18]. The eigenstates can be labelled by 
the total number of excitations N = X^( a ] a j + a j+ a j-) 
which is a good quantum number with [TV, H t ] = 0. In 
Fig. [2l the spectra for the low-lying eigenstates at the 
strong-coupling limit g ^> J and A = are shown. At 



FIG. 1: (Color online) Schematic circuit for the resonator 
array. Inside the dashed box: resonators where photons are 
generated; outside the box: interface to outside circuits. 

J = without the tunneling, the eigenstates are prod- 
uct states of the polariton modes in a single JC Model 
fl8| . At finite (but weak) J, the tunneling Hamiltonian 
is treated as a perturbation. For TV = 0, we have the 
ground state |^o) = IIj|0 For N — 1, the lowest 

four states \ipi lT n) with 1 < m < 4 are superpositions 
of the four lower polariton states L^-J^o). The tun- 
neling Hamiltonian lifts the degeneracy of these states 
with energy splittings hJ. Other states for N = 1 are 
the higher polariton modes T^-J^o) which are separated 
from the lowest states by an energy ~ hg. For N = 2, 
the lowest six states 1^2,™) for 1 < m < 6 are superpo- 
sitions of the states l)-^l) k -J^o) with j ^ k and contain 
two excitations occupying different sites. The tunneling 
Hamiltonian generates energy splittings of ^/2hJ among 
these states. All other states for N = 2 either involve 
excitations in the higher polariton modes or involve two 
excitations in the same site, such as 1^2,™) = ^\ 2^0) fc> r 
7 < m < 10, with an energy ~ hg above the lowest six 
states. 

Entangled photon pairs can be generated in the JC Lat- 
tice by spectroscopic techniques in the strong-coupling 
limit. The effective nonlinear interaction together with 
the tunneling J lifts the degeneracy of the eigenstates 
and makes it possible to generate novel quantum states 
by applying microwave pulses with selected phases and 
frequencies. Assume a monochromatic driving — 
hY^=i( e j e ~ lUJdta ] + ^je lUJdt aj) being applied to the res- 
onators with the amplitudes ej and frequency uod- This 
pulse can generate transitions between eigenstates differ- 
ing by one excitation number. When \ej\ <C J, g, only 
resonant transitions are important. Furthermore, we can 
choose the relative phases and amplitudes of (ei, 62, 63, €4) 
to design the allowed transitions. Consider the eigenstate 
at A = 

|^2,3) = - 4i L li)l^>/^ (2) 

which contains two excitations in the lower polariton 
modes. This state is an EPR state with spatially- 
entangled polaritons if we map the modes Li^ and L^^ 
to one polarization and the modes L 2 ,i and L^^ to the 



FIG. 2: (Color online) (a) Eigenenergies at J = A = 0. The 
degeneracies are labelled to the right, (b) Selected eigenstates 
at finite J and A = 0. Solid arrows are the designed tran- 
sitions to generate 1^2 ,3) and dashed arrows are off-resonant 
transitions. 

orthogonal polarization [12). Below we construct a two- 
pulse process to generate this state with high fidelity. 
Initially, with Huq ^> fc^T in the cryogenic environ- 
ment, the system is prepared to the ground state |^o) by 
thermalization. The first pulse has uniform amplitudes 
(ei, 62, 63, 64) = e(l, 1, 1, 1) with e <C J and the pumping 
frequency aid = ojo — g + J . This pulse induces a resonant 
transition and hence a Rabi oscillation between |^o) and 
the single-polariton eigenstate ^1,4) = (S^J,i'0o)/2, as 
is indicated by the solid arrow in Fig. EJb). After apply- 
ing the pulse for a duration of ir/2\f2e, the state is then 
pumped to | ^1,4) . Although nonzero transition elements 
exist between the state |^i,4) and some two-excitation 
states such as ^2,2), these transitions are suppressed 
by the off-resonances ~ hJ when e <C J. The second 
pulse has opposite phases in the odd and the even sites 
with (ei, 62, 63, 64) = e(l, — 1, 1, — 1) and the pumping fre- 
quency uod = uq — g — J. Among the states 1^2,™) 
for 1 < m < 6, the only nonzero transition element 
induced by this pulse is (^2,3 I^I^m) = No cou- 
pling is induced by this pulse between the states ^1,4) 
and l^o) ? and hence there is no transition back to the 
ground state. Although nonzero transition elements ex- 
ist between ^1,4) and other N = 2 states with higher 
energy such as ^2,7) and between ^2,3) and some of the 
N = ?> states, these transitions are again suppressed by 
the off-resonances ~ hg. This pulse then induces a Rabi 
oscillation between ^1,4) and 1^2,3) • After being applied 
for a duration of 7r/2e, this pulse generates the entangled 
polariton state 1^2,3)- 

The entangled polaritons can be transformed into entan- 
gled photons by adiabatically switching the JC Lattice 
to the "superfluid" regime. To achieve this, we slowly in- 
crease the detuning to enter the dispersive regime with 
A ^> g. During this process, the lower polariton modes 
adiabatically evolve into photon Fock states and the non- 
linear interaction is effectively turned off. The entan- 
gled state ^2,3), as an eigenstate, adiabatically evolves 
into (a\a\ — ajaj) |^o)/\/2 which contains two spatially- 
entangled photons. The adiabatic condition for reach- 
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ing this state requires that the increase of the qubit 
energy satisfies duo z /dt <C 4g 2 , and hence the duration 
of the adiabatic process satisfies St > (A/4g 2 ) pjj. 
For our system, this condition can be readily met when 
1/J < St < 1/e. Meanwhile, this state can also be writ- 
ten in the momentum space. In the "superfluid" regime, 
the single-particle eigenmodes are 
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with energies uoq — 2J for mode C2, co>o for modes c\ and C3, 
and uo + 2 J for mode C4. It can be shown that the above 
entangled state can be written as (c\c\ — c\c\)iI)q / \/2 and 
is hence entangled in the momentum space as well. By 
constructing coherent interface with outside circuits, the 
entangled photons can be distributed to outside modes 
and used in quantum information protocols. 

The success of the proposed scheme relies on three crite- 
ria: 1. the energy scales are well separated with e <C J, g\ 
2. the decoherence rates of the resonators ft and qubits 
7q satisfy ft, j q <C e; 3. the fabrication errors of the sys- 
tem parameters can be compensated. For superconduct- 
ing resonators |4|, |l9j , strong coupling to qubits reaching 
g/2ir = 100 MHz has been demonstrated. The tunneling 
J can be obtained as capacitive coupling between adja- 
cent resonators. For a superconducting coplanar waveg- 
uide resonator, the voltage at position x is 



Vj(x) = + aj) cos(2nx/L) (4) 

with L being the length of the quasi-one-dimensional res- 
onator. With a coupling capacitance C\ and resonator 
capacitance C r , we derive the coupling as [20] 

Ci^(L)VS-+i(0) = fuv ^-(a]a J+1 + a] +1 aj) (5) 

which yields hJ = huo(Ci/C r ). The tunneling can be 
adjusted in a large range by varying C\. For example, 
ujo/2tt = 10 GHz, d = 10 fF, and C = 2pF, we have 
J/2tt = 50 MHz. 

For finite driving amplitude e, the off-resonant transitions 
occur with the probability ~ \he/SE\ 2 depending on the 
off-resonance SE. For example, as is shown in Fig. [2j 
the second pulse in this scheme induces transition be- 
tween 1^1,4) and ^2,7) with the probability ~ |e/0.6g| 2 . 
Decreasing the driving amplitude can decrease the off- 
resonant transitions. However, it also means longer pulse 
duration and larger influence from the decoherence of 
the resonators and qubits. The damping rate of the res- 
onators can be as low as ft/27r ~ 10 kHz given a quality 
factor Q = 10 6 from recent experiments. The decoher- 
ence rate of the qubits can be less than 7 g /27r ~ 100 kHz. 



At A = 0, the decoherence rate of the polariton modes 
is 7 P ~ (ft + 7 g )/2 which can be less than 50 kHz. In 
the dispersive regime, the decoherence rate is mostly de- 
termined by the damping rate of the resonator. A driv- 
ing amplitude of e/2ir = 5 MHz can hence be adopted to 
achieve high fidelity in the generation of entangled states. 

To test this analysis, we numerically simulate this scheme 
using a master equation approach and calculate the fi- 
delity of the final density matrix: F = (^2,3 \pf ^2,3) 
with the target state ^2,3) • For simplification, we use 
the energy spectrum in Fig.[2jb) in our simulation, which 
can be treated as a Bose-Hubbard Model for the lower- 
polariton modes with damping rate j p . The effective 
self-interaction is defined as U = (E2- —2Ei_) /hin terms 
of the polariton energies [2] and U = (2- V2)g at A = 
when the entangled photons are generated. The other en- 
ergy levels in Fig. [2](a) are omitted. This model, though 
simplified, captures the dominant effects that reduce the 
fidelity: the off-resonant transitions and damping. The 
damping terms in the master equation are the standard 
Lindblad form for modes Lj : i [21]. In Fig. [3](a), the fi- 
delity increases rapidly with U in the regime U < J. 
For U > J, the fidelity stops increasing with U and 
reaches a saturation value that does not depend strongly 
on the driving amplitude. This agrees with our previous 
analysis that the fidelity is limited by the off-resonance 
SE ~ min(ftj, hU) when j p <C e. In Fig. [3^b), we find 
that the fidelity increases with the driving amplitude e at 
first, but soon reaches a maximum and starts to decrease 
as e further increases. The initial increase of the fidelity 
is due to the shortening of the pulse durations and hence 
the lessening of the damping as e increases. When e fur- 
ther increases pass the maximum, the off-resonant tran- 
sitions gain a larger probability which eventually makes 
the fidelity decrease. Our result shows that the fidelity 
can exceed 0.9 at 7 P - 100 kHz (j p /J - 2 x 10" 3 ). 

The parameters of the superconducting devices can 
spread by a few percent due to fabrication errors. The 
magnitude of the deviations in the resonator and qubit 
frequencies can be comparable with the coupling and 
the tunneling constants. Mostly coming from the in- 
accuracies in the Josephson energy, the deviations in 
qubit frequency can be avoided because the qubits are 
made of SQUID loops instead of single junctions. The 
qubit frequency as well as the detuning A is hence tun- 
able. Below we will show that by adjusting the detuning, 
the effect of the deviations in the resonator frequency 
and the deviations in the coupling constant on the en- 
ergy of the Lt j^o) modes can also be fully compen- 
sated. The energy of the L^ -J^o) modes is E\-/h = 
ujo + A/2 — y // A 2 /4 + g 2 [2] . Assuming the deviations Scuo 
in the resonator frequency and Sg in coupling strength, 
we find that the detuning needs to be 

A=[(5g + g) 2 -(5uj +g) 2 }/(8u; +g) (6) 
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FIG. 3: (a) Fidelity versus U/J at e/J = 0.02, 0.04, 0.1 from 
top to bottom and "y P /J = 2 x 10 -4 . (b) Fidelity versus e/J 
at 7 P /J = 0, 2 x 10" 4 , 2 x 10" 3 , 0.01 from top to bottom 
at U/J = 2. (c) Solid (dashed) curves: detuning (effective 
U/J) versus frequency errors Sujq at coupling errors Sg — 
O.lg, 0, —O.lg from top to bottom (from bottom to top). And 
J/2tt = 50 MHz. 



in order to recover the energy E\-/h = ujo — g for the 
lower polariton mode at A = 0. This detuning is plotted 
in Fig.[3](c) for various Suoq and Sg. Note that the energy 
of the T^-J^o) modes can not be compensated simulta- 
neously. As a result, the interaction U after the com- 
pensation varies with Suo and Sg and can be decreased 
for a negative Sojq. To ensure a large U to protect the 
fidelity of the pumping scheme, we can choose the lowest 
resonator frequency in the array as the base frequency 
ujo so that Sujo > for all other resonators. 

In summary, we presented a scheme to generate entan- 
gle photons in coupled superconducting resonators. Our 
scheme explores the Kerr nonlinearity in a JC Lattice, 
the circuit symmetry, and the spectroscopic techniques 
to pump entangled photon pairs with high fidelity. We 
numerically studied the robustness of this approach in 
the presence of decoherence. The idea can also be easily 
generalized to semiconductor system. When applied to 
multiple resonators, this idea can shed new light in the 
quantum engineering of novel states of microwave pho- 
tons (m. 
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